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Abstract 

We propose definitions of homogeneity and projective equivalence for 
systems of ordinary differential equations of order greater than two, 
which allow us to generalize the concept of a spray (for systems of 
order two) . We show that the Euler-Lagrange fields of parametric La- 
grangians of order greater than one which are regular (in a natural 
sense that we define) form a projective equivalence class of homoge- 
neous systems. We show further that the geodesies, or base integral 
curves, of projectively equivalent homogeneous differential equation 
fields are the same apart from orientation-preserving reparametriza- 
tion; that is, homogeneous differential equation fields determine sys- 
tems of paths. 
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1 Introduction 



This paper is concerned with systems of ordinary differential equations 

■ ■ ■ • • d'^v^ 

of order n + 1, n > 1, in m dependent variables , i = 1,2, . . . ,m and one 
independent variable x. Second-order systems of this kind {n = 1) have 
been studied intensively using methods of differential geometry, by taking 
advantage of the fact that the equations are those satisfied by the integral 
curves of a vector field of a certain type, a so-called second-order differential 
equation field (often abbreviated to SODE) on the tangent bundle TM of 
an m-dimensional differential manifold M, or some open submanifold of 
it. Higher-order systems may be studied by analoguous methods, where 
now the differential equation vector field in question lives (in the case of 
equations of order n-|- 1) on the bundle T^M of n- velocities, in other words 
n-jets at zero of curves in M (whose domains contain zero), or some open 
submanifold of it. This is the approach adopted here. (For a recent account 
of the geometric theory of differential equations, covering both second- and 
higher-order systems, and containing many references, see [1]. However, this 
reference is mainly concerned with aspects of the theory not covered here, 
namely generalized connections and related matters.) 

In the theory of second-order differential equation fields an important 
role is played by a class of fields which are homogeneous in a certain sense. 
Let ToM be the slit tangent bundle of M {TM with the zero section deleted). 
Let A be the Liouvillc vector field, that is, the infinitesimal generator of di- 
lations of the fibres of TqM. A second-order differential equation field T 
on ToM is called a spray if it satisfies [A,r] = F. The canonical spray of 
a Finsler space, whose integral curves project onto geodesies with constant 
speed parametrization, is a helpful example. If F is a spray then the distri- 
bution S)"*" spanned by F and A is involutive. Any second-order differential 
equation field F' contained in 2)+ is said to be projectively equivalent to F, 
and its base integral curves are obtained from those of F by reparametriza- 
tion. Slightly more generally, for the purposes of this paper we shall say 
that a second-order differential equation field F is homogeneous if the dis- 
tribution 2)+ spanned by F and A is involutive. It is not difficult to show 
that if F is homogeneous in this sense then there is a projectively equivalent 
second-order differential equation field which is a spray. In Finsler geom- 
etry the canonical spray is the Euler-Lagrange field of the energy, which 
is a regular Lagrangian. But instead of the energy one may consider the 
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Finsler function itself as a Lagrangian. Because of its assumed homogeneity 
it is parametric: that is to say, the corresponding action integral is invariant 
under reparametrizations. Such a parametric Lagrangian determines, not a 
single Euler-Lagrange field, but a projective equivalence class of second-order 
differential equation fields, which in the Finslerian case is the projective class 
of the canonical spray. 

Our purpose in the present paper is to generalize the concept of homo- 
geneity so that it applies to differential equation fields of higher order. For 
this purpose we introduce the jet group of order n, L", which is defined as 
follows. Consider the local diffcomorphisms (j) of M, defined in a neighbour- 
hood of zero and satisfying (/)(0) = 0. Then is the set of n-jets at zero 
of such local diffcomorphisms; it is a group under composition, and even a 
Lie group. Moreover, it acts to the right on T^M. Let 2) be the distri- 
bution spanned by the generators of the action. Wc say that a differential 
equation field P is homogeneous if the distribution 5)+ spanned by F and 5) 
is involutive. This generalizes the concept described above for second-order 
differential equation fields because when n = 1, 25 is just the one-dimensional 
distribution spanned by A. (A preliminary attempt to examine homogene- 
ity of higher-order differential equation fields is to be found in [2]; however, 
this paper is again concerend mainly with connection theory.) 

Section 3 below is devoted to the definition of homogeneity and its im- 
mediate consequences. 

With this definition wc show that the other properties of homgcneous 
second-order differential equation fields discussed above also generalize. The 
extremals of a Lagrangian of order n satisfy differential equations of order 
2n. If the Lagrangian is parametric, so that it satisfies the so-called Zermelo 
conditions, then its extremals are determined only up to reparamctrization. 
We show in Section 4 that, provided it is sufficiently regular, a parametric 
Lagrangian determines a projective equivalence class of homogeneous differ- 
ential equation fields. This result gives us the reassurence that our definition 
is not vacuous. 

In Section 5 we show that the base integral curves of projectively equiva- 
lent homogeneous differential equation fields differ only in parametrization. 

Section 6 of the paper contains some examples, while Section 2 is devoted 
to the requisite geometry of jet manifolds and jet groups. 

We shall use the Einstein summation convention for coordinate indices 
such as i\ other sums will be indicated explicitly. 



3 



2 Geometrical background 



We consider differentiable manifolds of class C°° which are Hausdorff , second- 
countable and (unless otherwise specified) connected. 

Let M be such a manifold with dim M = m and with local coordinates 

(y*). Consider curves 7 : (a, b) — )• M, where G (a, b), and let T^M be the 
set of n-jets at zero {jq t}) with coordinates (j/*) for < r < n given by 



It is a standard result that T^'-M is a manifold with dimT^M = (n + l)m, 
and that the maps r„ : T"M M and t„ „/ : r"M ^ T^'M (n' < n) given 

by 

rnOo"7)=7(0), V„'Oo7)=io'7 

are fibre bundles. 

The manifold T^+^M has a natural identification with a submanifold of 
TT"'M, obtained because the map t ^"7 defines a curve 

■■ i ^io(7 0Tt) 

in T^M, where : M — )• M is the translation x t-^ x + t. We therefore obtain 
the inclusion 

given in coordinates by y^^i 1-^ y*. We may think of this inclusion as a sec- 
tion of the pull-back vector bundle t*_^_^^^TT"-M T^'+'^M, in other words 
as a vector field along the projection r^+i^n! with such an interpretation it 
is the total derivative 

n Q 

The projection r„+i,n : T'^'^^M — > T^M is in fact an affine bundle, modelled 
on the vector bundle VTn,n-i ~^ T'^M of very vertical tangent vectors on 
r"M. The affine action is just addition in the fibres of TT^M T'^M: 

{yU0,...,0,z') + {yUyi,...,ytr-i,ytr) ^ {yUyl---,yi-i,z' + yi) 
gives rise, using the identification yl = on T'^^^M, to 

(.yUz') + {yUyii+i) ^ iyU^" + yti+i)- 
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We shall make use later of the vertical endomorphism on T^M. This is 
the type (1, 1) tensor field S given in coordinates by 

" d 

it is canonically defined, and generalizes the well-known vertical endomor- 
phism or tangent structure on TM. It is not difficult to show that for any 
1-form a on r"M, 

S{dTa) - dT{Sa) = r*+i_„a 

(where in the first instance S is the vertical endomorphism on T'"-^^M, in the 
second the one on T^^M). This may conveniently be written Sd'Y — d^S = 1, 
with the pull-back map understood (note that we assume action on 1-forms; 
there is a more complicated formula for action on forms of higher degree) . 

A section f : T"M T'^+^M of the affine bundle t„+i,„ is caUed a 
differential equation field; in coordinates it is indeed a differential equation 

y^+i = ^'(y^,•••,y;)• 
A curve 7 : (o, 6) — ^ M is a geodesic of the equation if J""'"^7 = T o J"7. 
Composing V with the inclusion T'^'^^M TT^M gives a vector field V on 
T'^M of the particular form 

n— 1 n rs 

every integral curve 7 : (a, h) — t- T^M of T is of the form 7 = j^j for some 
geodesic 7 of T; that is to say, the geodesies of T are the base integral curves 
of r. We use the term 'differential equation field' to refer to F as well as T. 

Now consider local diffeomorphisms ^ of M defined in a neighbourhood 
of zero and satisfying (f)(0) = 0. Let L" be the set of n-jets at zero {Jq of 
these local diffeomorphisms; this is a group under composition, 

Jo<^i •io'/'2 = Jo(<^io<^2) 

and is a Lie group. As a manifold it has two connected components. The 
component of the identity is a subgroup L""*" of index 2 containing jets of 
local diffeomorphisms satisfying (?!>'(0) > 0. The map 

j>0^(<^'(0),<^"(0),...,<^(-)(0)) gM- 



5 



is a global coordinate system on (and L""*"). We can obtain an explicit 
formula for the product in these coordinates by using an expression for 
the nth derivative of a composition of functions. This is Faa di Bruno's 
formula [4], which we take in the form involving the Bell polynomials: 

(eor/)H(x) = f^e^'^)(r?(a;))S; {rj' {xW [x), . . . ,ri^^+'-^\x)) ; 

r=l 

here ^ and -q are functions of some variable t with (in our case) ^(0) = 77(0) = 
0, and B^{rii,ri2, ■ ■ ■ , ?/n+i-r) is a polynomial in the n + 1 — r variables r)p, 
for which the following explicit formula is known: 



Bn{ril,V2,---,rin+l-r) 



n! /^2\'?2 / Vn+l-r ^''"+^" 



- 



where the sum is over all non-negative integers qi,q2, ■ ■ ■ , qn+i-7- such that 
qi + q2 + ■■■ + (In+i-r = r and qi + 2q2 + . . . + {n + 1 - r)qn+i-r = n. (This 
formula may easily be derived by talking ^(x) = x'^ and expressing ry(x) as a 
formal power series in x: 

ri(^x) = —x + —x + —x +•••; 

then {r\/n\)Bl^{r]i, r]2, ■ ■ ■ , 77n+i-r) is the coefficient of x" in the formal power 
series for {rj{x)Y . In principle could depend on rjp with p > n + l — r, but 
in practice, since ri{x) contains no constant term there can be no contribution 
to a;" in {r]{x)Y coming from the term rjpX^ /p\ if p > n + 1 — r.) To express 
multiplication in in terms of the global coordinate system introduced 
above, a vectorial notation is convenient. We denote by ^, 77 elements of 
M", considered as row vectors, and by 6(77) the (upper triangular) matrix 
whose p, q entry is Bq{r]i,r]2, ■ ■ ■)■ Then the product in L" is given by 

|77 = |B(r7); 

we must of course assume that .^1 and r/i are non-zero. Note that the ma- 
trix which determines multiplication in with n' < n is just the n' x n' 
submatrix of B(t7) in the upper left corner. For instance, taking n = 4 gives 



B(r?) 



r/f 3r/i?72 4771773 -|- 377^ 
ryf 677^772 
V r]f ) 
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There are obvious projections A„,i : ^ L?- and A^^ : L""*" 
given by they are homomorphisms because 

An,l(io '/'I • io -^s) = K,l{jQ{4>l o 02)) = jo{(t)l ° (t>2) = ioVl " jo(f>2- 

The kernel of A„^i is therefore a normal subgroup ii'" <l L". By inspection 
of the matrix B above in the case n = 4 we see that and are abelian 
but K'^ is not; indeed is non-abelian whenever n > 4. 

We may, in addition, identify with a subgroup of by mapping 
Jq(/) to jQ^(f,i{Q)-, where //^ : 1^ ^ is the multiplication diffeomorphism 
ljLs{t) = st for s 7^ 0. Thus L" may be regarded as a semidirect product 

XI ii'", using the action of on by conjugation. 

Furthermore, we may write = L'^ y\ J"', where J"" is the kernel of the 
homomorphism A„,2 '■ ^ given by j^cp i-^ j^cp. For the inclusion of 

in L" we take, as a representative of jg^, not the quadratic polynomial 
4){x) = ax + \hx'^ but the Mobius transformation 

4>{x) 



1 — a ^bx 



because the composition of two of these transformations is another map of 
the same form. By calculating the derivatives of (f) one finds a polynomial 
representative of j^cf), 

r=l 

(see [5]). Note that the kernels of the homomorphisms A^^^ do not give 
semidirect product decompositions when r > 2. 

We next obtain a basis for the Lie algebra T of the jet group in terms 
of the coordinates defined above. For this purpose we need the following 
lemma. 

Lemma 1. For given n, p = 1,2, . . . ,n, and r = 1,2, . . . ,n + 1 — p, 



dr)r 

unless p = n + 1 — r , when 



= 
(i,o,-,o) 



dB, 



drjr 



ni 

(i,o,...,o)~ ("-0!^!' 
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Proof. We consider the cases r = 1, r > 2 separately. We have 



dm 



{qi - l)!g2! • • • gn+i-p'l! 
,91-1 



X 



Vn+l-r 



{n + 1- r)\ 



gn+l-r 



where the sum is over ah non-negative integers qi,q2, . . . , 9„+i-r with > 
such that qi+q2 + - ■ .+qn+i-r = P and qi+2q2 + . . . + {n + l-r)qn+i-r = n. 
With ryi = 1, ?7r = for r > 2 we get a non-zero contribution to the sum 
only if = for r > 2, when we must have qi = p and qi = n. That is, 
there is no non-zero term in the sum unless p = n, and so 



dm 



(l,0,-,0) 







nl 



[ (n-l)!l! 



p ^ n 
p = n 



For r > 2 on the other hand 



dVr 



qi\- ■ ■ {qr - ly. ■ ■ ■ g„+i_p!r! 



X 



{n + l- r)\ 



In+l-r 



where the sum is now over all non-negative integers qi,q2, ■ ■ ■ , qn+i-r with 
qr>0 such that qi+q2+- ■ .+qn+i-r = P and qi+2q2+. . .+{n+l-r)qn+i-r = 
n. With r/i = 1, ?7s = for s > 2 we get a non-zero contribution to the sum 
only if (7,. = 1, qg = for s > 2, s ^ r. We must therefore have qi + 1 = p 
and + r = n, that \s, p = n + 1 — r and qi = n — r. Then 



dBf, 



(l,0,-,0) 







n! 



{n — r)\r\ 



p ^ n + 1 — r 
p = n + 1 — r 



□ 



We denote by (J'", r = 1, 2, . . . , n, the left-invariant vector field on L" 
which takes the value 

r\ — 
' dyr 

at the identity. Here (yi,y2, • • • iVn) are the coordinates on L". The factor 
r! is included for later convenience. We now obtain an explicit expression 
for (5^ in terms of coordinates. 
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Proposition 1. 

Proof. Let y = {yi,y2, ■ ■ ■ , Vn) be a generic point of L". The identity e has 
coordinates (1, 0, . . . , 0). We have 

= -^yl*e<^e 

where Ly is left multiplication by y: Lyt] = y'B{r]). Then 

d 



^ s=l p=l 



(l,0,-,0) 



s 



n , 

S\ 



Es: d 
7 rr-^ys+i-r-K—- □ 
(s — rj!r! oys 



s=r 



Corollary 1. 



(r — s)5^~^^ ^ ifr + s<n + l 
otherwise 



□ 



One reason for introducing the factor r! in the definition of is to simplify 
the expression for the bracket. 

We can now make some observations about the Lie algebra T which 
correspond to properties of the Lie group L" mentioned earlier. In the first 
place, {6^,6^} span a subalgebra. Furthermore, for any p = 1,2, 
{6^ : r > p) is an ideal in P (since for any r > p and s > 1, r + s — 1 > p). 
Let Ip C P be the ideal just defined (so that in particular P = Ij*). Then 

for p>i, r/[;j ~ p-i. 

The ideals for p > 1 are nilpotent. It is known that for a nilpotent Lie 
algebra, exponentiation is a surjective map onto the simply connected Lie 
group of which it is the algebra; it follows in particular that exponentiation 
maps [2 onto i^" (see [5]). But it is possible to prove this directly, as we 
now show. The Lie algebra P consists of linear vector fields on M", and so 
exponentiation in the group coincides with matrix exponentiation. In the 
case of an element k = Ylr=2 ^^5^ °f '2 matrix K„ to be exponentiated 
is strictly lower triangular: its elements are given by 



' T <s 



(Kn)r 



I {s-l) 



kr+i-s r>s 
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Notice that the (n — 1) x (n 
the matrix K„_i correspond 
example, with n = 5 we have 



1) submatrix in the upper left corner is just 
the matrix K„_i corresponding to the element Ylr=2 ^rS^ of Ig"^. As an 



/ 



3! 









\ 







V §h f|fc4 §h §k2 J 










Hi., ill, 4[u 



Proposition 2. For any y = (1,7/2,^3, ■ ■ ■ ,yn)^ S K"- there is k G I2 *^c/t 
t/iat y = exp n. 

Proof. As a column vector in M", expK G L" is just (expK„)e, where e = 
(1,0,..., 0)^; that is to say, exp k is the left column of exp K„. It is easy to 
see that the first n — 1 entries in the left column of exp K„ comprise the left 
column of the (n — 1) x (n — 1) matrix exp K„_i, while the lower left corner 
of exp K„ is of the form n\kn +p{k2, k^, . . . , kn-i) where p is a polynomial in 
the indicated variables. Since K„ is strictly lower triangular, expK„ has Is 
down the diagonal. We have to show that for any y2) 2/3, • • • , 2/n we can choose 
k2,ks, . . . ,kn such that the left column of exp K„ is (1, 7/2, 2/3, • • • , Vn)'^- We 
proceed by induction. Thus we assume that we can find k2,k-^, . . . ,kn-i 
such that the left column of expK„_i is {I,y2,y3, ■ ■ ■ ,yn-i)^- With these 
values of ^2, ks, . . . , kn-i we set kn = {yn — p{k2, k^, . . . , fc„_i))/n!: then the 
left column of exp K„ is (1, y2, 2/3, • • • , Z/n)^ as required. We have 

^"P^^=(k2 

SO the first step certainly works. □ 

We next turn to the relationship between and certain vector fields on 
M. An element 6 of [" determines a one-parameter subgroup of L", namely 
t exp{tS): observe first that the integral curve of S through y G is 
t y ■ exp(t(5). Now let be a one-parameter group of diffeomorphisms 
of M such that 0f(O) = 0. Then cpt = Jq '^t is evidently a one-parameter 
subgroup of L" whose infinitesimal generator, a vector field x on L", is left 
invariant and therefore an element of P. For any y G L^, Xy is the tangent 
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vector at t = to the curve t y ■ But the sth component of y ■ is 
given by 



p=l 











\ dx 


o' 


5 • • 




■' dx^+P-^ 



Note that since (f)o{x) = x, the argument of at t = is (1, 0, . . . , 0), while 
d fd^(l)t\ fd(t)t\ (FX 

(0,0) \ / (0,0) 



dt V dx^ 



dx"- 



where X{x)d/dx is the infinitesimal generator of (j)t (a vector field on R). 
On differentiating the formula for (y • ^t)s with respect to t at t = we 
obtain 



X 



In fact 



dPX 



(s — dxP 



d dPX 



6P. 



^ ^ ^1 dPX 
^ p\ dxP 



p=i 



SP 



is a linear map from the space of vector fields on M vanishing at the ori- 
gin onto P, such that ^{Xd/dx) depends only on j^X. In particular, 
^{x'^d/dx) = (5^ for r = 1, 2, . . . , n while ^{x^d/dx) = for r > n. Now 



d 



d 



dx^ dx 



(s — r)x 



r+s-i ^ 
dx^ 



so $ is an anti-homomorphism. 

Let p be the Lie algebra of vector fields on R whose coefficients are formal 
power series in j;, and let p" be the subalgebra of those vector fields which 
vanish to order n at 0, that is, whose coefficient begins with Then 
p° is the subalgebra of formal power series vector fields which vanish at 0; 
and for n > 0, p" is an ideal in p°. If we think of $ as a map p*' — [" it 
is a surjective anti-homomorphism with kernel p", and therefore defines an 
anti-isomorphism p^'/p" I". To put things another way, we can realise 
P as the space of vector fields on R whose coefficients are polynomials of 
order n which vanish at the origin, with bracket the negative of the ordinary 
bracket of vector fields followed by truncation at order n. This alternative 
realisation is found elsewhere in the literature. The fact that the bracket is 
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related to the negative of the ordinary bracket is not too surprising when 
one recalls that if the diffcomorphism group of R is regarded as an infinite 
dimensional Lie group, its Lie algebra is the space of vector fields on M with 
compact support, but with bracket the negative of the ordinary vector field 
bracket. 

The group L" has a right action on T"M given by composition of 
jets, 

a„ : L« X T"M ^ T^M, (io>,io7) ^ Jo (7 o 0); 

the action is fibred over the identity on M. It restricts to the 'regular' 
submanifold T^M (that is, of ra-jets Jq 7 where 7 is a curve with 7'(0) 7^ 0, so 
that 7 is an immersion near zero) because composing a diffeomorphism with 
an immersion gives another immersion. We may use Faa di Bruno's formula 
recursively to see that the restricted action is free, using the coordinates 
ivl), ^ <i <n, defined on complete fibres over M: by regularity at least one 
coordinate y\ must be non-zero at any given point, so if anijo4'j JqI) = Jo 7 
we sec successively that i^'(O) = 1 and then that 4>"{0) = 0, 0'"(O) = 0, 

so that = The orbit space of T^M under the action 

has a manifold structure (and is, indeed, a Hausdorff manifold); we shall 
denote it by PT'^M. The subgroup L"+ acts in the same way, and its orbit 
space will be denoted by P+T^M; this is a double cover of PT'^M. Let 
Pn : T^M PT^'M and p+ : T^M P+T^-M be the projections. 

In the present work we shall in effect be interested in the circumstances 
when a differential equation field on T^M (more accurately, on T^M) 'passes 
to the quotient' to determine a line-elcmcnt field on PT^M or an oriented 
line-element field on P^T'^M. In view of the identification of T^~^^M with 
an affine sub-bundle of TT^M over T^^M, it is of some interest to consider 
the identification of PT"+^ with a submanifold of PTiPT^-M). One may do 
this, for instance, by using the action of the tangent group TL" on TT^M 
to give T{PT'^M), and then the action of on the open submanifold 
To{PT"'M) to give PTiPT^-M). In fact it is known that pr"+i is an affine 
bundle over PT'M, even though PT{PT''M) -)■ PT'M certainly does 
not have an affine structure, as its fibres are compact (they are projective 
spaces): see, for instance, the discussion in [6], which uses the fact that the 
kernel of the homomorphism L"+^ — )■ L" is abelian. Similar arguments hold 
in the oriented case. 

We shall need to know the fundamental vector fields on T"'M of the 
action of L". We denote by A** the fundamental vector field corresponding 
to € P . Then A** is the infinitesimal generator of the one-parameter group 
i2exp(t(5'-) (where R denotes the right action). In terms of the coordinates (y*) 
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on T'^M, the action of ?7 G L" is given by 

r 

p=i 

By a by now familiar type of argument invoking Lemma 1 we obtain 
Proposition 3. 

In particular A^, which corresponds to the infinitesimal generator of dila- 
tions of R, is given by 



Corollary 2. 

[A^A 



•SI 



(r-s)A''+^-i ifr + s<n + l 
otherwise 



Corollary 3. In terms of the vertical endomorphism S 

1. a vector field T on T^M is a differential equation field if and only if 
S{T) = A\- 

2. A'^+i = SiA-") = 5^(Ai) = 5^+i(r). 

Wc shall not use the results of Corollary 3 directly, but wc mention 
them because they provide an alternative starting point for the geometrical 
analysis of differential equation fields which may be found elsewhere in the 
literature. 

We shall however use in Section 4 a rcsTilt which is related to item 1 
of Corollary 3, which we now explain. We denote by i the operator of 
interior product of a vector field with a form, so that for example for a 
1-form a, ixa = a{X). We extend this notation to apply to the total 
derivative: if a is a 1-form on T^M then ixct is a function on T^~^^M. Then 
iT(<S'Q;) = ('^n+i,n'^)' which may conveniently be written i^S = i^i, with 
the pull-back map understood. Note that we assert only that this formula 
holds for action on 1-forms (there is a more complicated formula for action 
on forms of higher degree). Similarly, in relation to item 2, we have (again 
for action on 1-forms) iTS^ = i/^r . 
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3 Homogeneous higher-order systems 

We now consider differential equation (vector) fields on T^M. 

Definition 1. A differential equation field T on T^M is homogeneous if the 
distribution spanned by T and the A*" is involutive. 

Proposition 4. The differential equation field T is homogeneous if and only 
if there are functions A**, r = 1, 2, . . . , n such that 

[A\r] =r + AiA", [A^^] =rA'-i + A''A", r = 2,3,...,n. 

When they exist, such functions must satisfy the following consistency con- 
ditions, where 1 < r, s < n, r ^ s: 

A\y) - A'-iX^) = (n + 1 - r)A''; 
A"(A^) - A^(A") = (r - s)X'+'-\ r + s<n + l; 
A'^(A^) - A^(A'') = -(n + l)(r - s), r + s = n + 2- 
A'^(A^)-A"(A'^) = 0, r + s>n + 2. 

Proof. A straightforward coordinate calculation shows that [A^,r] — F and 
[A'',r] — rA''"^ are very vertical. Thus in order for S)+ to be involutive 
there must be functions A'' such that 

[ASr] =r + AiA", [A",r] =rA''-i + A'^A", r = 2,3,...,n. 

The consistency conditions follow from the Jacobi identities. □ 

In terms of coordinates the conditions for 

r = Vy;+i-- + r-— 

^ + dyl dyl 

to be homogeneous are 

A\r) = {n + l)r + nW'yl A^V) = J^lill^y^^^-. + nW^yi 

If r is a differential equation field such that S)"*" is involutive, so is r+/xA" 
for any function ji on T^M, and of course it belongs to the same involutive 
distribution. 
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One cannot, in general, demand that for n > 3 a differential equation 
field r satisfies the strong conditions 

[A\r]=r, [A^^] = rA'-^ r = 2,3,... ,n; 

the Jacobi identities would be inconsistent with the bracket relations for the 
A^'. (This point is discussed more fully in [2], where examples of equations 
which do satisfy the conditions for n = 2 are given.) Recall that {A^, A^} 
span a subalgebra of 2); it is easily verified that the conditions 

[A\r] = r, [A2,r] = 2Ai 

are consistent. We shall now show that for any homogeneous T it is possible 
to find n such that if F = T + //A" then F satisfies the conditions above 
(as well as [A'',r] = rA"-^ + A'' A" for r > 2); in other words, if F is 
homogeneous, without loss of generality we may assume that = = 0. 
We need a couple of lemmas. 

Lemma 2. There is a submanifold S of T^M of codimension 1 such that 
A^ is transverse to S and A^ is tangent to it. There is a submanifold S' of 
S of codimension 1 such that A?'\s is transverse to S' . 

Proof. Let g be any Riemannian metric on M. Then gijUiul is a well- 
defined function on T^M (because the y\ transform as the components of 

tangent vectors), and gijy\y{ > on T^M. Let (p{y) = \J gijy\y{, and 
S = {y e T^M : ip{y) = 1}. Then A^{(p) = and A'^{ip) = 0, so A^ and 
A^ are respectively transverse and tangent to <S. Next, let (p' = ^T{(p) for 
any differential equation field F on T^M. We have 

A\^') = ^A2(F(<^)) = iF(A2(<^)) + A\<p) = ^, 

since [A^,F] — 2A^ is very vertical. Now 

^vv' = h^yWivl + 9ijyiyi = myliyi + ^idvi)^ 

where the F^,- are the connection coefficients of the Levi-Civita connection 
of the metric from which it is clear that S' = {y ^ S : <p'{y) = 0} is a 
codimension 1 submanifold of 5, and A^l^ is transverse to it. □ 

Lemma 3. Let X he a complete vector field on a manifold M , with 1- 
parameter group 4>t, and S a codimension 1 submanifold of M transverse 
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to X such that {^t(<S) : t E M.} = M. Let f be any smooth function on 
M, k any constant, and zq any smooth function on S. Then there is a 
unique smooth function z on M such that X(z) + kz = f and z\s = zq. In 
particular, if f = and zq = then z = 0. 

Proof. One can integrate the differential equation 

dz 

- + kz = fiM^)), 

where x is any fixed point of <S, by the integrating factor method. □ 

(The restriction that k should be constant covers the situation encountered 
below, but clearly more general equations could be considered.) 

Theorem 1. LetT be a homogeneous higher-order differential equation field 
on T^M , so that S"*" is involutive. Then there is a higher-order differential 
equation field F = F + //A" G 2)+ such that 

[A\T]=T, [A2,r]=2A\ 

Proof. Let 5 be a codimension 1 submanifold of T"M such that A-*^ is trans- 
verse to S and A^ is tangent to it, and S' a codimension 1 submanifold of 
<S such that A^l^ is transverse to S'. 
We may write 

[A\f] =f + A^A", [A2,f] = 2Ai + A2A'*. 

Then with F = f + ^A", 

[A\T] = F + {A\n) + (2 - n)ii + A1)A". 

By Lemma 3 wc can solve the equation A^(/x) + (2 — n)/x = — A-^, assigning 
the value of /x on 5 arbitrarily. Also, 

[A2, F] = 2Ai + (A2(/x) + P)A". 

Now A^ is tangent to S, and jj, is undetermined on that submanifold, so we 
may solve the equation A^(//) = — A^ there, assigning the value of /x on <S' 
arbitrarily. With satisfying these two equations, for F we have A^ = 
everywhere on T"Af and A^ = on 5. But from the consistency condition 
for A^ and A^ applied to F we have A^{X'^) = (n - l)X^ on T^M. It follows 
from Lemma 3 and the fact that A^ = on ^2 that A^ = everywhere. □ 
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We suggest, somewhat tentatively, that homogeneous higher-order dif- 
ferential equation fields such that [A-'^,r] = T and [A^,r] = 2A^ should be 
considered as generalized sprays. 

Finally in this section we observe that a homogeneous differential equa- 
tion field on T"M does indeed pass to the quotient under the action of 

Theorem 2. Let T be a homogeneous differential equation field on T^M. 
Then the involutive distribution D'^ defines an oriented line-element field 
on P+T'^M. 

Proof. The distribution 2)"'" evidently projects to a line-element field on 
p+j"i^M. To orient it we take, at any point Pnijol)^ the positive multiples 
of Pn*^j^'y foi' any representative point 7 G T'^'-M. □ 



4 Parametric Lagrangians 

To show that the theory described above is not vacuous, we devote this 
section to an important source of examples, the higher-order differential 
equation fields obtained from parametric variational problems. (Related 
results, but for second-order variational problems only, may be found in [8]. 
Multiple-integral parametric variational problems are discussed at length 
in [3], and indeed the methods used here are related to those used in that 
paper, but of course specialized to the case of a single integral.) 

Let L be a Lagrangian function, defined on T^M. The Hilbert form of 
L is the 1-form 

^ = y ±-^dPsp+'dL 

defined on T^"'~^M. This may be used to construct the Euler-Lagrange 
form of the Lagrangian, 

e = dL — d'Y'd 
defined on T^^M. By construction, therefore. 



n-l 



e = dL-dTY -^-^(E,SP+^dL = y L-^dF^SPdL. 
^ (p + l)\ ^ „ p! 

We may see that e is horizontal over M, because S'c/t = cfxS' -|- 1 (modulo 
pullback), so that 

S($^ = -|- pS^^ 



V 
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and therefore 

Thus Se is a collapsing sum, and indeed 



Sd^SP = dPSP+^ + pdF-^SP. 



s£ = y LJ^d^sp+^dL + y EL^dP^-^sPdL = ^-^d^^s^+^di = o 

p=0 ^ p=l ^ 

as L has order n. 

Lemma 4. The Euler- Lagrange form e vanishes along the extremals of any 
fixed- endpoint variational problem defined by L. 

Proof. An extremal of such a problem is a curve 7 in M such that 

rb 



j irircxnLdt = o 

J a 



for any vector field X on M vanishing at the (fixed) image of the endpoints 
of [a,b], where X"' denotes the prolongation of X to T^M. Then 

^\rindix^L)dt- 



j {rjr£xnLdt= f {r^r{ix-dL)dt+ j 

J a J a J a 



but 

rb 



f {rir{dixr^L)dt= f d{{r^rixnL)dt=[{rjrixnL]i = o 

J a J a 

as X vanishes at the endpoints, so that 
f\riTCxr^Ldt= f\riT{ixr^dL)dt= [\f^jnix2nie + dT^))dt. 



But now, 

rb 



[ {f''^nix^ndT^)dt= [ (f"7)*(f^T^x2n-l^?)c^^ 

J a J a 

= j\{{f^-'l)%ix^n-.d))dt=[{f--'^)%ix^n-.d)]\=Q 

J a 

because vector field prolongations commute with total derivatives, and the 
pull-back of cZt to M is d. Thus 

''\f''intxe)dt 



j {rircx^Ldt= [ {f^^r{ix2^e)dt= [ 

•J a J a J a 



because e is horizontal over M. As X is arbitrary, it follows that (J^"7)*e = 
0. □ 
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Now suppose that L is positively homogeneous; that is, that it satisfies 
the Zermelo conditions 

A\L) = L, A'^(L) = (r>2). 

Such a Lagrangian is called a parametric Lagrangian, and the geodesies 
of the corresponding variational problem are invariant under orientation- 
preserving reparametrizations [3]; in other words, they are paths. We shall 
show that (subject to suitable regularity conditions) they are also the geodesies 
of homogeneous differential equation fields. In order to investigate this, we 
first establish a technical lemma. 

Lemma 5. The Hilhert form 'd and the Euler- Lagrange form e satisfy 

= L, i^di^ = —e 
(omitting the pull-back maps). 

Proof. From iq^dT = (^t^t and iTS^ = i^p we obtain 

iTd^SP+'^dL = d^iTSP+^dL = dP.i^p+idL 
J i/^i dL = L ip = 0) 



so that 



In addition, 



(p > 0) 



n-l 

= y / ' iTd^SP+'dL = L. 



dxi? = diT^ + ztc?'!? = dL + i^d-d 
so that, from the definition of e, 

£ = dL- d-rd = -i-rd-d. □ 

Proposition 5. The characteristic distribution of the 2-form d'd satisfies 

{r,A\A2,...,A2"-i} ckerdi? 

whenever T is a differential equation field on T^M whose geodesies are ex- 
tremals of the variational problem with fixed endpoints defined by L. 
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Proof. We remark first that i? satisfies = and Ca^i^ = for 1 < r < 
2n — 1 ([3], Proposition 6.1 and Theorem 6.4). It fohows that 

lArd'd = C^^-d — diA'r'& = 0, 

and so {A\ A^, . . . , A^"-!} c kcvdi}. 

Now let r be a differential equation field on T^^~^M, so that 

f : T'^'^-^M T^^'M, 

and let T : T^^'-^M i{T^''M) C TT^'^-'^M be the corresponding vector 
field. 

Suppose in particular that the integral curves of T (which must, neces- 
sarily, be prolongations J^"^^7 of curves 7 in M) are such that 7 is always 
an extremal of the variational problem. Each point of T^^~^M must lie on 
such an integral curve, and we may suppose (by translation in the domain) 
that the point in question is J^"~^7(0), in other words that the point may 
be written as Jq"^ 7- 

The integral curve property means that r(jQ"~"^7) = jo"7, so that 

r..„-l^ = (iof)(j2"-l^)=i(j-2n^); 

thus for any vector ^ G T.2n-i^T^'^~^M we have 

= 

because 7 is an extremal of L so that £j-2n^ = 0. This calculation also uses the 
fact that the contraction is just the inclusion map i : T^'^M TT^'^~^M 
in disguise. □ 

We say that the Lagrangian is regular if the characteristic distribution 
of d-d has dimension 2n. Any differential equation field in this distribution is 
an Euler-Lagrange field of the variational problem. Since the characteristic 
distribution of a closed 2-form is involutive we have the following theorem. 

Theorem 3. Any Euler-Lagrange field of a regular parametric Lagrangian 
is homogeneous. 
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5 Projective equivalence and reparametrization 



We now return to the study of general higher-order differential equation 
fields. 

Definition 2. Let T, T' be homogeneous. Then if T' — T = //A" for some 
function n we say that V and V are projectively equivalent. 

Projective equivalence is an equivalence relation. The projective equiva- 
lence class of a homogeneous differential equation field consists of all the dif- 
ferential equation fields in the involutive distribution 2)+. We may summa- 
rize the result of Theorem 1 by saying that every projective equivalence class 
of a homogeneous differential equation field contains a generalized spray. 

We have seen that the Euler-Lagrange fields of a regular parametric La- 
grangian consist of a projective equivalence class of homogeneous differential 
equation fields. So we should expect in this case that the geodesies of projec- 
tively equivalent homogeneous differential equation fields should define the 
same paths, that is, differ only in parametrization. In this section we shall 
show that this is true for any projective equivalence class of homogeneous 
differential equation fields, whether or not they come from a parametric 
Lagrangian. 

In fact we shall prove two results about projective equivalence in homo- 
geneous systems: firstly, if V is projectively equivalent to F then a geodesic 
of r' is a reparametrization of a geodesic of F; secondly, the jet group acts 
on differential equation fields in such a way as to map any homogeneous 
field to one which is projectively equivalent to it. 

Theorem 4. Let V he a homogeneous differential equation field, and V 
another which is projectively equivalent to it. Letj' be a geodesic ofV. Then 
there is a geodesic 'jofT such that 7' is obtained from 7 by an orientation- 
preserving reparametrization. 

Proof. By assumption the distribution is involutive, and thus integrable. 
Let £ be a leaf of S)"*". Then F' is tangent to £, and any integral curve 
of F' which meets S, lies in it. If 7' is a geodesic of F' then f-j' is an 
integral curve of F': let £ be in fact the leaf of 2)+ containing (the image 
of) J" 7'. Since !l) C 2)^, £ is invariant under the action of L""*", the identity 
component of the jet group L", and consists of the orbits of points on j^^' 
under the action of Under projection r„ : T"M M, we have of 

course t„ o ^"7' = 7', and for any g G L"+, Tn o g = Tn. Thus r„(£) is the 
1-dimensional submanifold of M which is the oriented path of 7'. Now F 
also belongs to S)"*", and so likewise any integral curve of F which meets £ 
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lies in it. Take any point of £, and the integral curve of T through that 
point: it is of the form j^'j where 7 is a geodesic of T. Then 7 = t„ o J^j 
lies in r„(£), that is to say, 7' and 7 have the same oriented path. Since 
we are restricted to (yl) 7^ 0, we must assume that 7 and 7' have nowhere 
vanishing tangent vectors. It follows that one is a reparametrization of the 
other, by a reparametrization that preserves orientation. □ 

Corollary 4. Let T be a homogeneous differential equation field of order 
n + 1. Let 7 he the geodesic ofT with initial conditions {y^,y\, . . . ,yn). Let 
{z\, . . . , z^) be the image of {y\, . . . , y^) under some element of L"^^ (so that 
{z\) is a positive scalar multiple of {y\))- Then the geodesic ofV with initial 
conditions (y*, . . . , zjj) is obtained from 7 by an orientation-preserving 
reparametrization. 

The second result is rather more complicated to prove, and requires some 
preliminaries. 

We shall need to distinguish notationally between the fundamental vector 
fields of the actions of on T"M and L"+^ on T^+^M. Let us denote by 
the vector field operating on T^M: that is to say 

With a bit of licence we can write, for r = 1, 2, . . . , n, 

n + 1)! , d 



while 



{n + l-rr^^-^-^dyl^, 
d 



In particular. 

Now the homogeneity conditions for an (n + l)th order differential equation 
field r, when expressed in terms of the F*, are 

Ai(r) = (n + i)r + n\\\\, A;(r) = ^^'l\^}\^^ U2-r + n\yy\. 

Comparison with the formulae for ^n+iiVn+i) suggestive. 
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Consider the section F of T^+^M T^M corresponding to T, given by 
Vn+i = r*(yo' yi> • • • 1 Vn)- We have, on imf, with r = 2,...,n, 

K+iiVn+i - n = (n + iWn^, - r) - nlX'yl = -nlX'yl 
A;+i(y;+i-r) = -n!AVi 

One way of restating this is that for r = 1, 2, . . . , n the vector fields 



r 



\ 

(strictly speaJiing one should write r*^^ ^A'' rather than just A^) are tangent 
to imF. Set 

\r 

The A'' are well-defined vector fields on r"+^M. 
Proposition 6. For r, s = 1,2, . . . ,n 



[A^A^] 

Proof. We have 



(r - s) A''+*-^ i/ r + s < n + 1 
otherwise 



[IS ,is \- i (nTT) 

while 

FA^ A*l — TA'' A* 1 I i'^ni^ ) ~ ^n(^ )) An+1 

(the A*" are functions on T"M). The stated results follow easily, using the 
consistency conditions on the A'". The one case which is not entirely straight- 
forward isr-|-s = n-|-2, when we have 



[A^ A^] = A^+i] + ^-^iS^lJ-—J^A^_,, 



(A;( A-)-A- (AO),,+i 
(n + l) 



since in this case A^(A*) - A^(A'') = -{n + l)(r - s). □ 
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That is to say, when T is homogeneous there is a representation of the 
algebra I" on imF. 

Notice also that for r > 2, [A;^+}, A*"] = 0, since [A^^+^.A^+^j = and 

A;:i}(A'-) = o. 

Let (j) be an orientation-preserving local diffeomorphism of R leaving 
fixed. Let (j)n be the diffeomorphism of T'^M given by the action of Jq G 
L". We have r„_|_i^„ o (f)n+i = 4'n Tn+i,n- Let F be an (n + l)st-order 
differential equation field, T : T^M T^^^M the corresponding section. 
Since r„+i,n o f = idT^M-, 'Pn+i ° ^ o is also a section of rji+i^„, and so 
defines a new (n + l)st-order differential equation field, say P^. Suppose 
that r is any homogeneous [n + l)st-order differential equation field. We 
shall show that for any (j), there is a function V(f, on T^M such that 

V\{yi)=V\yi) + u^{yi)y\. 

Note that when this holds V^ = V + v^fyA^, so Fff, is projectively equivalent 
to r. From the point of view of representation by sections one might say 
that f(t, is obtained by sliding f along the integral curves of A^^^J. To be 
more precise, recalling that T^+^M is an affine bundle over T"'M modelled 
on the bundle of very vertical tangent vectors on T^M, is the translate 
of F by the very vertical vector field f^A^ljlJ. 

Theorem 5. Let T be a homogeneous differential equation field and let (j) 
he an orientation-preserving local diffeomorphism of M leaving fixed. Then 
F^ is projectively equivalent to V. 

Proof. Note first that ^(t>i°<t>2 — (J^<l>i)(j>2- So if for some 4>i,(f)2 and every 
homogeneous F, it is the case that F,^-^ is projectively equivalent to F and 
F02 is projectively equivalent to F then F<^jo(^2 projectively equivalent to 
F. We pointed out in Section 2 that the finite order jet group L'^ is the 
semi-direct product of and K"^, the subgroup consisting of the jets of 
diffeomorphisms with derivative at equal to the identity, and further that 
every element of the latter is the exponential of some element of its Lie 
algebra. So the result about F^, or equivalently F^, will be proved if it can 
be proved when = kx {k a positive constant) and when (j)n+i is the 
exponential of an element of (A5^_|_]^) for r = 2, . . . , n. 

Consider first the case where (j){x) = kx. Then (j)n{y\.) = (k^yl). We 
have Ajj(F*) = (n + 1)F* -|- nl\^y\, which integrates to give 

g-(n+i)tp(er-tyi) _ p(y^) = ri\ f e-^'\\e''yi)y\ds. 

Jo 
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Then 

where ^ ^ 

Jo 

Next, let nt be the 1-parameter group generated by ^^^=2 ^p^n+i ^o'^ some 
constants kp. Now 

n+l n I ^ n \ 

p=2 p=2 \ ^ p=2 ) 

Let be the flow generated by X]p=2(^p^^)' ^'^'^ ^^^^ generated by 
TJp=2(kp^n), SO that 

l^t ° Tn+l,n = Tn+l,n O ht = T„+l,„ O Kj. 

For any (y*) G T"M let T{t,yl.) be the solution of the first-order ordinary 
difi^erential equation 

such that r(0, y*) = 0. We claim that 

where (j)^^^ is the 1-parameter group generated by A^^}. Consider the 
(dim T'^M -|- l)-dimensional submanifold S of rj*+^M consisting of points 
0^+i(f(y^)), t G M, (y^) G r„"M. Since Ep=2(^P^^) is tangent to imf, 
Kt(f(y;)) G imf for all t, and so (?!)^+^^i^(Kt(f(y*)) G 5 for all t. We can 

use coordinates (t, y*) on 5; with respect to these coordinates 0"'''^ is just 
translation of the first coordinate, and 

^;!+i.)(«t(f(y;)) = (r(i,y;),Kt(y;)). 
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Moreover, since [A";]l];,A^] = the coordinate representation of co- 
incides with its coordinate representation on imF, which is the same as 
the coordinate representation of AJ^. Thus the tangent vector to the curve 
1 1-)- ^"^^j^(/cj(f (y*)) at t, in coordinate form, is 

d " 

p=2 

But d/dt is the coordinate representation of A^^J on S, and we conclude 
that t 0"J^i^(Kt(f(y*)) is an integral curve of the vector field 

( n \ n n+1 

^+%=2 / p=2 p=2 

Since (/)^+i^,^(Ko(f(y;)) = f(y^), we see that t ^ cP^^+^^.^ktir {,/,)) is the 
integral curve of ^^^=2 ^p^n+i starting at T{yl), that is, it coincides with 
Kj(r(y*)) as claimed. 
It follows that 

But fit maps imf to itself, and Tn+i,n okt = Rto Tn+i,n, so Kt(f = 
t{yi). Thus 

Kt{r{R_M))) = €tt!^.,(yi))(^(yr)) = {y'rX{yl) + r{t,R^M)y\). 

This establishes the required result with Kt = 4>n+i ^(piVr) = ^i^^ 4>n{yl))- 

□ 

This result shows that a projective equivalence class is invariant under 
the action of the group L"+ on T^M. 



6 Examples 

Take M to be Euclidean space of dimension m and the to be Euclidean 

coordinates. With the scalar product etc. denoted in the usual way, consider 
the third-order differential equation field defined by 

pi/ s ( 'i\yi?\y2? + {yi ■ y2? \ i .^iyi^m) i 
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This is easily seen to be homogeneous, with Ai = A2 = 0. 

We now explain the geometrical significance of this third-order system. 
Denote the differential equation field by F. It is easy to see that 

r(|yi|2) = 2(yi-y2) 

2(yi • y2) 



r(yi • y2) 



It follows that r is tangent to the submanifold on which \yi\ = 1, yi ■ 2/2 = 0. 
This consists of the 2-jets of curves parametrized by Euclidean arc-length. 
The restriction of T* to this submanifold is just — |y2pyi- The corresponding 
vector differential equation, considering the as the components of a vec- 
tor r with respect to an orthonormal frame and using overdots to indicate 
differentiation with respect to arc- length s, may be written 



Along any solution curve of this differential equation the 2-plane spanned 
by r and r is constant, and so the curve is a plane curve. Moreover, is 
the curvature of the curve, and 



d_ 
ds 



rP = 2r • f = -2|rPr • r = 0. 



So the curve is a plane curve of constant curvature, that is, a circle (or if 
the curvature is zero, a straight line). 

Now the submanifold \yi\ = 1, yi ■ y2 = in this case is just the subman- 
ifold S' of Lemma 2. It follows that every point of T^M can be written in 
the form <i)ti{4't2('^)) ^'^^ some z e S', where 4>l and (f)f are the 1-parameter 
groups generated by A-*^ and A^. So in particular every point of T^M IS 
the image by some jet map of a point of S' . It follows from Corollary 4 
that the homogeneous system we started with also has the property that its 
geodesies are circles. 

There will of course be other homogeneous differential equation fields 
with the same property, namely those projectively equivalent to T. For 
example we could take the simpler system with 

r (yi,y2) = 3 j/2- 

\Vi\ 

The corresponding field is homogeneous, now with = (yi • y2)/|yiP (but 
A^ = still). However, this field is not tangent to «S', so while it is true (by 
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Theorem 4) that its geodesies are circles, it is not so straightforward to see 
this. 

We next exhibit a fourth-order system in Euchdcan space derived from 
a parametric second-order Lagrangian. Consider the Lagrangian function L 
on T^M given (in terms of Euchdean coordinates (y*) as before) by 

L{y,yi,y2) = p-rg ■ 

mr 

This satisfies A^{L) = L and A^(L) = 0, so is parametric. In fact L is 
closely related to the (first) curvature of a curve, which may be considered as 
a function k on T^M: we have L = (see for example [9]). By deriving 

the Euler-Lagrange equations for L one obtains the projective equivalence 
class of fourth-order equations, necessarily homogeneous, containing the one 
given by 

j.i _ o(y2jJ/3) i , /^ 5|yip|y2p -35(j/i ■y2f + 8\yi\'^{yi-y3) \ ^ (yij/2) i 

ml V 2|yi|4 J 

For this particular representative F we have = 0, but is nonzero, and 
so of course is A'^. We find that F(|yi|, F(yi • 1/2) and F(yi • 7/3 -|- |y2p) are 
contained in the ring generated by these functions, which means that F is 
tangent to the submanifold S" of T^M where \yi\ = 1, yi-y2 = ?/i-y3+|y2p = 
0, whose points consist of jets of curves parametrized by arc-length. On S" 
we obtain the differential equation 

r'= — 3(r • r)r — ||rpr 

(in terms of arc- length parameter). This may be more succinctly written as 

^(•r + i|r|^r)=0. 

This is the Euclidean version of an equation discussed recently by Matsyuk 
in the context of 'Zitterbewegung' [7] (eq (38) with ^ = and R = 0). 
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